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Abstract. This Minireview summarizes the recent efforts to solve for- 
ward and inverse problems as they occur in different branches of funda- 
mental and applied magnetohydrodynamics. As for the forward prob- 
lem, the main focus is on the numerical treatment of induction pro- 
cesses, including self-excitation of magnetic fields in non-spherical do- 
mains and/or under the influence of non-homogeneous material param- 
eters. As an important application of the developed numerical schemes, 
the functioning of the von-Karman-sodium (VKS) dynamo experiment 
is shown to depend crucially on the presence of soft-iron impellers. 
As for the inverse problem, the main focus is on the mathematical 
background and some first practical applications of the Contactless 
Inductive Flow Tomography (CIFT), in which fiow induced magnetic 
field perturbations are utilized for the reconstruction of the velocity 
field. The promises of CIFT for fiow field monitoring in the continu- 
ous casting of steel are substantiated by results obtained at a test rig 
with a low melting liquid metal. While CIFT is presently restricted to 
fiows with low magnetic Reynolds numbers, some selected problems of 
non-linear inverse dynamo theory, with possible application to geo- and 
astrophysics, are also discussed. 



1 Introduction 

When a moving electrically conducting fluid comes under the influence of a magnetic 
field, an electromotive force is produced that drives electrical currents in the fluid 
virhich, in turn, modify the original magnetic fleld. This general principle applies both 
to the case that a magnetic field is applied from outside, as well as to the case that it is 
exclusively produced inside the fluid. In the latter case, we speak about self-excitation, 
or homogeneous dynamo action. It can only occur if the magnetic Reynolds number 
Rm = fiaLV, i.e. the product of magnetic permeability fi, electrical conductivity 
(T, size L and velocity scale V of the moving fluid, is signiflcantly greater than one. 
This condition is well fulfllled in the convecting fluid layers of cosmic bodies, such 
as the Earth [T] or the Sun [2, but it is much harder to achieve in the liquid metal 
laboratory. Yet, the threshold of self-excitation has been reached in some dedicated 
liquid sodium experiments in Riga [3^, Karlsruhe and Cadarache [51, and it is 
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targeted in a few more places around the world [B]. If, however, Rm is not large 
enough for self-excitation, a magnetic field applied from outside will still be modified 
by the moving fluid. This effect can be utilized for inferring the velocity field from 
the externally measured induced magnetic field [7]. 

In this Minireview, we will deal with both effects: induction under the influence of 
applied magnetic fields, and self-excitation. We will discuss two mathematical formu- 
lations, the first one, the differential equation approach, being based on the induction 
equation, the second one, the integral equation approach, relying on Biot-Savart's 
law. In either case, particular attention will be paid to the correct treatment of the 
boundary conditions for the magnetic field, which is a notorious and non-trivial prob- 
lem in non-spherical geometries. We will exemplify both methods by treating forward 
problems related to the French von-Karman-Sodium dynamo experiment for which 
we will show that the use of magnetic propeller materials plays a decisive role for the 
functioning and the mode selection of this type of dynamo |8I9) . 

The integral equation method is then applied to various inverse induction problems 
in 2D and 3D. We will show that the Contactless Inductive Flow Tomography (GIFT), 
as we call it now )7i, has a promising potential for the online flow monitoring in 
technologies such as steel casting or crystal growth. The much harder, since non- 
linear problem of inverse dynamo theory is only touched upon by discussing some 
simplified models. 



2 Induction and self-excitation: Mathematical basics 

In this section we will discuss two different approaches to induction and self-excitation, 
the first on based on differential equations, the second one on integral equations. 



2.1 Differential equation approach 

2.1.1 The induction equation 



The evolution equation for the magnetic field B in a fluid of electrical conductivity 
a and relative permeability fir can be derived from Ampere's law, Faraday's law, the 
divergence- free condition for the magnetic field, and Ohm's law in moving conductors: 

Vx^=j (1) 

VxE^-- (2) 

V • B = (3) 
j cr(E + v X B). (4) 

In Eq. (2) we have skipped the displacement current because in good electrical conduc- 
tors the quasi-stationary approximation can be applied. We will also restrict ourselves 
to the case without external currents. That way, taking the curl of Eq. (1) and Eq. 
(4), and inserting Eq. (2), we obtain the induction equation for the magnetic field: 



xB+ ^'' xB — -VxB . (5) 

firfiQO- fir fJrf^OO- J 

This is the general form in which spatially varying conductivities and relative 
permeabilities are taken into account jB], as it will indeed be considered in the later 
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treatment of the VKS dynamo experiment. If the material parameters are constant 
in a given volume, the induction equation reduces to the better known simple form: 

^ = V X (v X B) + AB . (6) 

The right hand side of Ec^. (jB]) describes the competition between the diffusion and the 
advection of the magnetic field. Comparing the diffusion time-scale with the kinematic 
time-scale gives the most important dimensionless number of magnetohydrodynamics, 
the magnetic Reynolds number Rm := pL^HraLV 

Depending on the actual flow pattern, the values of the critical Rm^ at which 
the field starts to grow, are usually in the range of 10^... 10^ p]. Most flows in cosmic 
bodies in which Rm is large enough will act as dynamos, although there are a number 
of anti-dynamo theorems excluding too simple structures of the velocity field or the 
self-excited magnetic field [lOj . 

Throughout this paper, the velocity field will be considered as given or, in the 
inverse problems, as an unknown which we will try to reconstruct. That means, in 
turn, that no attempts will be made to solve the Navier-Stokes equation in the MHD 
regime. 



2.1.2 Simulations in spheres 

The easiest geometry for treating the induction equation, together with the demand 
that B is curl free in the exterior and vanishes at infinity, is the spherical geometry. 
Fortunately, planets and stars are in first approximation spherical so that dynamo 
problems for those bodies are comparably easy to solve. 

Considering the solenoidal character of magnetic fields, V • B = 0, and the incom- 
pressibility of liquid metals, V ■ v = 0, both the flow and the magnetic field can be 
represented in a toroidal-poloidal decomposition: 

v(r) = V X V X s(r)f + V x i(r)r, (7) 
B(r) = V X V X S'(r)r + V X T(r)r. (8) 

Following the BuUard-Gellman formalism [11] the scalar functions s,S,t,T can then 
be expanded in terms of the spherical harmonics Ym.ii'd, (p): 

5(r,z9,(^) =^5,,™(r)rz,™(i?,^) and T(r, z9, (^) = ^ Tj,™(r)l^,„(i9, (^). (9) 

where we denote the degree and order of the spherical harmonics by I and to, respec- 
tively 

By inserting Eqs. (|7I8I) in Eq. ©, multiplying with the conjugate complex f), 
integrating over Mdip, and applying the orthogonality relation of spherical harmon- 
ics, one can derive a coupled system of equations for the defining scalars Si^mir) and 
Ti^m{r) which are second order in r. More details of this procedure can be found, e.g., 
in |12) . What is interesting in our present context is that the insulating boundary 
conditions, which require V x B = in the exterior, lead to two boundary conditions 
for the scalar fields Si^m and Ti_„i at the outer radius R that are diagonal in I and to: 

dSi.^jr ^R) 

7{ y ISi.rnir ^ Rj = 

or 

Ti,m{r = R) = 0. (10) 
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It is this diagonality that makes dynamo simulations in spherical geometry quite 
simple, and the spherical harmonics decomposition was indeed the method of choice 
for many simulation of cosmic dynamos, but also for the optimization of spherical 
dynamo experiments such as the MDX experiment in Madison |13] . 

2.1.3 Treating the boundary conditions for non-spherical domains 

The convenient treatment of the boundary conditions, which are diagonal in the 
degrees and orders of the spherical harmonics, does not extend from spherical to other 
geometries. For that reason the simulation of liquid sodium dynamo experiments, most 
of them working in cylindrical geometry, needs significantly higher numerical effort. 

For rather compact machines like the Karlsruhe dynamo experiment, with its 52 
spin generators and a cylindrical aspect ratio close to one, it appears reasonable to 
utilize the convenience of spherical harmonics expansion by embedding the actual 
cylindrical vessel into an enclosing virtual sphere and to assume the space between 
this enclosing sphere and the real cylinder to be filled with matter of low conductiv- 
ity. Choosing conductivity factors of 1, 100, and 1000 between the cylinder and the 
embedding, a convergence towards a realistic critical Rm was obtained in ^1^ . 

While this embedding method works reasonably well for cylinders with an as- 
pect ratio close to one, its application to very long dynamos is inconvenient due 
to the rather slow convergence. Hence, for the Riga dynamo another method has 
been employed [15116117] which is based on the fact that taking the limit /iu — s- 
the induction equation reduces to a Laplace equation. While in the dynamo region 
the induction equation is solved, the Laplace equation is solved in the exterior by 
a pseudo-relaxation method. The solutions in both domains are than matched by 
appropriate interface conditions for the magnetic field. The scheme has been used 
extensively for the prediction, optimization and data analysis of the Riga dynamo 
experiment and the resulting growth rates and frequencies were in remarkable agree- 
ment with the experimental ones. 

In comparison, the utilization of simplified vertical- field conditions only (i.e., van- 
ishing tangential fields at the boundary, B*^'^ — 0), turned out to underestimate the 
critical Rm of the dynamo by some 20 per cent, an error margin that is not acceptable 
for predictions of expensive liquid sodium experiments. 

Another method to treat insulating boundary conditions has been developed 
in [18119120] . The scheme uses a spectral finite element method and was applied 
successfully to several problems in connection with the VKS experiment [8121122] 
and to precession-driven dynamos in cylindrical domains [23] . The starting point is 
the fact that the divergence-free and curl-free magnetic field in the non-conducting 
exterior can be expressed by a gradient of a scalar magnetic potential, according 
to 'Bext = — V<?. Then, the Finite Element approximation for the magnetic field 
H = B/iJ-r in the conducting interior is matched to the solution of the outer Laplace 
equation by using an Interior Penalty Galerkin method. Typically, the outer spherical 
insulating domain is made ten times larger than the inner domain so that the concrete 
form of the boundary conditions at the very outer boundary becomes irrelevant. 

A possibility to avoid the solution of the Laplace equation in the exterior is to uti- 
lize the mathematically equivalent boundary element method (BEM). This procedure 
has been proposed in two papers by Iskakov et al. [24125] . and was recently modified 
and applied to various dynamo problems [26127] . We will present this method in some 
more detail, since it seems to be most efficient for the simulation of dynamo problems 
in arbitrary domains. 

The method starts with the normal component of the magnetic field at the bound- 
ary as obtained from a finite volume method in the interior. The unknown tangential 
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components of the magnetic field at timestep (n + 1) can then be derived by a ma- 
trix operation on a vector composed of the normal components of B at all surface 
elements of the computational domain. 

The base of the BEM is once again the possibility to express the magnetic field B 
as the gradient of a scalar magnetic potential ^ which fulfills the Laplace equation: 

B = -V<P with A<P = 0, <P^ 0{r-'^) for r ^ oo. (11) 

The computation of B on the boundary therefore requires the integration of A(I> = 
0. Consider a volume f2 that is bounded by the surface F where ^/on denotes the 
derivative in the normal direction: '9/9n = n • V so that 9n<? = — i?" yields the normal 
component of B on F. We can write Greens second identity for the scalar magnetic 
potential and a test-function G: 

JiGA<P-^AG)df2^J {Gf^-<^'-^)dF (12) 
n r 

Using Z\<? = 0, the potential is determined by the integral expression 

*W = /(G(r,r02|^-*,O^^2£^).r(,', (13) 

r 

where G(r, r') is the Green's function, or fundamental solution, given by 

and therefore fulfilling Z\G(r,r') — (5(r — r'). In principle, Eq. (IT^ is only valid for 
V ^ F . For all source points on the boundary the integration domain is enlarged by a 
small half sphere with the radius e and a short calculation of the limit e — > results 
in the boundary integral equation (see e.g. [26]): 

i.<,./(.,,.,2g_,,,2^).,,,, 

From (fTS]) the tangential components of the magnetic field on the boundary = 
e,- • B = — ■ V<?(r) are computed by: 

B- = 2jer- (^{v')Vr^^^^^ + S"(r')V,G(r, r')) dr(r'), (16) 
r 

where e-T represents the tangential unit vector on the surface element dF{r'). 

A discretization of the equation system ([T5|) and yields an algebraic system of 
equations which allows for the computation of the (unknown) tangential components 
of the magnetic field. The discretization of F in boundary elements is realized by the 
tessellation resulting from the finite volume discretization on the domain surface. 



2.2 Integral equation approach 



An alternative way to deal with the non-local boundary conditions of MHD is the 
integro-differential formulation developed by Meir and Schmidt |28I29I30] . Using the 
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current density as the primary electromagnetic variable, it is possible to avoid artificial 
boundary conditions and fully account for the electromagnetic interaction between 
flow region and surrounding. The point is that the magnetic field has then to be 
expressed in terms of the current by Biot-Savart's law. This method is presently 
being used in a realistic simulation of liquid metal experiment on the so-called Tayler 
instability, a current-driven kink type instability |31| . 

Interestingly, this method can be extended into an iterative scheme, by which 
not only induction effects, but also self-excitation can be treated. Starting with an 
applied field Bq, the authors of [52] search for solutions of the total magnetic field 
B = Bo -t- Bind = Bo -I- J^kLi '^k, with Bfc = 0{Rm'^). Actually, from a given Bfc 
they compute the electromotive force u x B^, from which the Poisson equation for the 
electric potential of the next iteration Aipk+i = V • (u x B^) is derived. From this, the 
electric current of the i + 1 iteration is obtained according to — —V'^k+i + vLX Bfc, 
and from this current the magnetic field of the i + \ iteration is obtained via Biot- 
Savart's law 



Despite the convergence problems of this method when approaching the critical 
Rm of a dynamo, it is quite helpful to illustrate the various induction processes that 
enter the dynamo mechanism. 

With view on Eq. p7)) it is tempting to ask why not to consider Biot-Savart's law 
as an integral equation for the magnetic field which has to be solved self-consistently. 
Assuming an infinite domain of homogeneous conductivity, such an approach has been 
utilized early by Gailitis |33) for the dynamo action of a pair of annular vortices and 
later by Dobler and Radler [34] for various dynamo problems. 

An extension of such an integral equation method which incorporates the effect of 
boundaries can already be found in the book of Roberts [33]. However, real use of the 
method was only made later in the papers |36l37j . With [38] this integral equation 
method was then extended to the case of unsteady magnetic fields which requires an 
additional equation for the vector potential A which is related to the magnetic field 
by B = V x A. The general formulation of the integral equation approach is the 
following: 



where Bo is the externally applied magnetic field (which might be zero) , b the induced 
magnetic field, u the velocity field, A the vector potential, (/) the electric potential 
and n denotes the outward directed unit vector at the boundary 5'. For a steady 
velocity field, the time dependence of all electromagnetic fields can be assumed to be 




(17) 




(20) 



(18) 



(19) 
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~ expXt. In (TTBl - Bni) the temporal variations of the velocity field are ignored, which 
is justified when the typical diffusive timescale of the magnetic field is shorter than 
the timescale of velocity variations (which usually is the case in the context of liquid 
metal experiments). 

We have to distinguish three different cases: For non-zero Bq, and below the self- 
excitation threshold, the imaginary part of A is simply the angular frequency of the 
applied and also of the induced magnetic field. For Bq = the equation system 
(|18H20p represents an eigenvalue equation for the unknown time constant A whose 
real part is the growth rate, and its imaginary part the angular frequency of the 
eigenmodes. For Bq = and A = 0, we need only Eqs. (US]) and which then 
reduce to an eigenvalue problem for the critical value of the velocity u at which the 
dynamo starts to work. 

Note that the integral equation system (jl8fl20l) is by far not the only possible one. 
There are indeed other possible schemes, one of them starting from the Helmholtz 
equation for the vector potential which leads, however, to a nonlinear eigenvalue 
problem in A (see |34j ) . while the above scheme is a linear eigenvalue problem in A 
which has many advantages when it comes to the numerical simulations. 

The numerical implementation of the integral equation approach was applied to 
various dynamo problems [38 .39 40.41) . and helped also to understand the under- 
performance of the first VKS experiment in terms of the detrimental role of (rotating) 
lid layers behind the impellers [39] . 



3 Forward problems: One example 

In this section we will invoke both the differential equation approach and the inte- 
gral equation approach for understanding the so-called von-Karman-dynamo (VKS) 
experiment in Cadarache [5], France. 

The working principle of the VKS dynamo is illustrated in Fig la. The flow of 
liquid sodium in a cylinder with diameter 412 mm and height 524 mm is produced by 
two counter-rotating impellers, each consisting of a massive disk and 8 curved blades. 
The resulting flow is of the so-called s2+t2 type, which comprises two poloidal eddies, 
indicated by "s2", with radial inflow {+) in the equatorial plane, and two counter- 
rotating toroidal eddies, indicated by "t2". This flow topology had been early shown 
to be a good candidate for dynamo action |^42 . In connection with the optimization 
of the VKS dynamo experiment, Marie, Normand and Daviaud j43j had studied a 
rather simple analytical test flow of this topological type, which is expressed in the 
following way: 

= r (1 - r)2(l + 2r)cos(7rz) (21) 

= 4er(l - r) sin(7rz/2) (22) 
-u^ = (l-r)(l + r-5r2)sin(7rz). (23) 

In order to study dynamo action for this MND-flow, with correctly implement- 
ing the vacuum boundary conditions, we had utilized both the differential equation 
approach (with solving the Laplace equation in the exterior) as well the integral equa- 
tion approach [35]. For the simplest case without any additional layers around the 
dynamo, and a toroidal/poloidal ratio e = 0.7259, the structure of the resulting mag- 
netic eigenfield is illustrated in Figs. lb,c. In Fig. lb the field lines of the equatorial 
dipole are seen, while the isosurfaces of the magnetic field energy (Fig. Ic) show the 
banana-like structure typical for this equatorial dipole eigenmode. 



8 Will be inserted by the editor 




(a) (b) (c) 

Fig. 1. The VKS dynamo, and its numerical simulation by the integral equation method, (a) 
Sketch of the experimental setup. The von-Karman swirling flow is driven by two counter- 
rotating impellers, (b) Magnetic field lines of the equatorial dipole field resulting for the 
simplifying MND flow, (c) Iso-surface of the magnetic field energy. 



A somewhat surprising results of [32] was that different kinds of layers around the 
flow field region can have opposite effects on the critical magnetic Reynolds number. 
While a (stationary) side layer works usually in favor of dynamo action, any layer 
(stationary or, even worse, rotating) behind the impellers impedes dynamo action 
significantly. This observation was one of the motivations for the VKS group to replace 
the non-magnetic impeller material by soft-iron, connected with the idea that such 
a magnetic material might shield the detrimental induction effects from behind the 
impellers. 

As a result of this modification, dynamo action was indeed observed [5], and 
subsequent experiments have brought about a remarkable variety of interesting phe- 
nomena, including field reversals, excursions, and bursts [44j. Such effects had been 
neither observed in the Riga nor in the Karlsruhe experiment and thus instigated a 
strong interest in the geophysics community, with a number of attempts to under- 
stand reversals of the geomagnetic field and reversals in the VKS dynamo on the 
common basis of low-dimensional models [45^ . 

Yet, the basic dynamo mechanism underlying such exciting results remained un- 
explained for some time. In contrast to all numerical simulations of the VKS dynamo, 
which had predicted a non-axisymmetric (m=l) eigenmode similar to that shown in 
Fig. 1, the actually observed eigenmode turned out to be more or less axisymmetric 
and dominated by the toroidal component. Furthermore, the experimentally observed 
critical Rm around 30 was significantly lower than the predicted one of about 50. The 
question arose, therefore, whether both effects could be explained by the utilization 
of magnetic impellers. 

In order to answer this question, we have utilized the combined Finite Volume/ 
Boundary Element Scheme exposed in 2.1.3. With this code we treated a slightly 
simplified model of the VKS dynamo, in which the bended thin blades were replaced 
by straight thick ones (see Fig. 2a). Our solutions show two different regimes that 
differ in the topology of the leading eigenmode. If the relative magnetic permeability 
of the disks exceeds a threshold of around 20, we observe an interchange of the leading 
eigenmode from the usual poloidal one to the toroidal one |5]. Yet, this magnetic disk, 
together with an assumed MND flow in the bulk of the fluid, does not lead to dynamo 
action, it just gives a toroidally dominated, but still decaying eigenmode (Fig. 2b). 
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(a) (b) (c) 

Fig. 2. Efltect of the magnetic impellers on the eigenmode selection in the VKS dynamo, 
(a) Assumed permeability distribution; isolevel at 87.5 percent of the peak value. In the 
fluid region fir is equal to unity. (b,c) Distribution of magnetic energy density and geometric 
structure of the magnetic field for fj,r = 20. (b) Rm — 30, a = 0: decaying axisymmetric 
mode with vanishing poloidal component, (c) Rm — 30, a — —1.5: Growing axisymmetric 
mode with poloidal component (after [9]). 

What is still needed for dynamo action is some amount of helical turbulence, param- 
eterized by the so-called a effect, which can well be assumed to exist in the space 
between the blades ^2 ■ This a effect is then able to close the dynamo loop by trans- 
forming toroidal field into poloidal field components [§]. The resulting eigenmode, 
which can now have a positive growth rate, is depicted in Fig. 2c. The experimen- 
tally observed eigenmode structure, with its strongly dominant toroidal component, 
is much in favor of our numerical simulation which explains the functioning of the 
VKS dynamo by means of strong suction of magnetic flux into domains with enhanced 
permeability (so called paramagnetic pumping |46| ) . together with some a effect con- 
centrated between the blades. This example shows that a very detailed treatment of 
induction effects, with a careful consideration of gradients of material parameters, 
might be essential for the understanding of experimental results. 



4 Inverse problems: Some examples 

Assuming a velocity field as given, and asking for the resulting induced (or self- 
excited) magnetic field, represents a typical forward problem. However, the direction 
of interest can also be inverted; assume the magnetic field to be known from mea- 
surements in the exterior of the fluid, what is then the velocity that produces this 
field? 



4.1 Inverse problems at low Rm 

Based on the integral equation system (I18ll20p in the time- independent limit (i.e. 
A ^> 0), we have developed the so-called Contactless Inductive Flow Tomography 
(GIFT) for the reconstruction of velocity fields from externally measured magnetic 
fields |38I48I49I50| . For the case of small Rm, which applies well to many industrial 
problems like steel casting or silicon crystal growth, the total field B = Bq -I- b under 
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the integrals in (jlSMlOp can be replaced by the applied magnetic field Bq alone. This 
replacement makes the inverse problem of reconstructing the velocity field from the 
induced electro-magnetic field b a linear one, which is much easier solvable than the 
full-fledged non-linear problem with B = Bq -I- b under the integrals. 

In a first step [48l49j . we had shown that, with a vertical field Sq.z being applied, 
the velocity structure of the fiow can be roughly reconstructed from the measurement 
of the induced magnetic field in the exterior and of the induced electric potential at 
the boundary of the fluid. There remains, however, a non-uniqueness concerning the 
exact radial distribution of the flow, a fact that can be made plausible by representing 
the fluid velocity by two defining scalars (for the poloidal and the toroidal velocity 
component), living both in the whole fiuid volume. Then it is clear that two quan- 
tities measured only on a two-dimensional covering of the fluid cannot give enough 
information for the reconstruction of the two desired three-dimensional quantities. 

Later, the reconstruction method was advanced to a completely contactless method, 
GIFT, to avoid the electric potential measurement at the fluid boundary [50]. The 
main idea was to apply the external magnetic field in two different directions (e.g. 
Bq z and Bq^x) and to utilize both corresponding sets of induced magnetic fields for 
the velocity field reconstruction. 




Fig. 3. Application of GIFT to a 3D flow driven by a propeller in a cylindrical vessel, (a) 
Schematic sketch of the experiment: M - Motor, C - Coils, P - Propeller, G - Guiding blades, 
H - Hall sensors. (b,c,d) Special case of upward pumping: (b) externally measured induced 
magnetic fields for applied Bo.z and (c) Bo,x, (d) velocity field as reconstructed from the (b) 
and (c). 

The goal of a first demonstration experiment [35] for GIFT was the reconstruction 
of a propeller-driven three-dimensional flow of the eutectic alloy GalnSn in a compact 
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cylindrical vessel with a ratio of height to diameter close to one (see Fig. 3a). In order 
to determine both the poloidal (in r and z) and the toroidal (in ip) flow components, we 
applied subsequently a vertical (-Bo,z) and a horizontal {Bq^x) magnetic field, either 
produced by Helmholtz-like coil pairs. In this experiment, the switching between 
the two fields occurs every 3 seconds, so that after 6 seconds all the magnetic field 
information is available for the velocity reconstruction. Note that this rather poor 
time resolution can be significantly enhanced, with a physical limitation given by the 
magnetic decay time, which is in the order of 0.05 s for the demonstration experiment. 
For each of the two applied fields, Bq z and i?o,a:: the induced fields are measured by 
Hall sensors at 48 positions, which are rather homogeneously distributed around the 
surface of the cylindrical vessel. The small ratio of 10^'^...10^^ between the induced 
and the applied fields demands for a very stable current source of the Helmholtz-like 
coils, for a very stable relative position of Hall sensors and coils, and for significant 
effort to compensate drift and sensitivity changes of the sensors. 

For the reconstruction of the velocity field we solve the Gauss normal equations 
that minimize the mean quadratic deviation of the measured from the modeled in- 
duced magnetic fields [35]. The intrinsic non- uniqueness problem concerning the de- 
tailed depth-dependence of the velocity is overcome by utilizing the so-called Tikhonov 
regularization by minimizing, in parallel to the magnetic field deviations, also the ki- 
netic energy of the flow or other appropriate quadratic functionals of the velocity. 
Mathematically, the Tikhonov regularization minimizes the total functional 

^[v] - Fbo.,, [v] + Fb,,^ [v] + F^,, [v] + Freg [v] (24) 

with 

Nb 



Fbo,. [-]=J2^ - b^"-\-]) ' (25) 

i=l ' 

Fbo,. [v] - E i i^as - bT'-'^ [v]) (26) 

. Nv 

div 

, Nv 

FregM^-^Y.^lAVk. (28) 

pen f.^^ 

The two functionals in (|25p and in (1^^ represent, for applied transverse field Bo,£c 
and axial field B0.2, respectively, the mean squared residual deviation of the measured 

(B ) (B ) 

induced magnetic fields h\ ^^^^ from the modeled fields h\ " [v] (that result from the 
solution of the corresponding forward problem) . Fdiv [v] (Eq. [27]) enforces the velocity 
field to be solenoidal, and Freg[v] (Eq. [^5)1 is the regularization functional which 
tries to minimize the kinetic energy. The parameters ai in (j25p and (|26p are the 
assumed a-priori errors for the measurement of the induced fields. The parameter 
adiv in (123) is chosen very small as it is a measure for the divergence the velocity 
solution is allowed to have. The parameter Upen in (|28p determines the trade-off 
between minimizing the mean squared residual deviation of the observed fields and 
minimizing the kinetic energy of the estimated velocity field. The normal equations, 
that follow from the minimization of the functional are solved by Cholesky 

decomposition. This numerical scheme can easily be extended to incorporate further 
a-priori information, e.g., the velocity components or mass flow rates known from 
other measuring techniques. 



12 



Will be inserted by the editor 



In the model experiment it was possible to distinguish clearly between upward and 
downward pumping of the propeller, with the rotational component being reduced 
by guiding blades in the case of upward pumping. For this case, Figs. 3b, c show the 
induced fields for applied i3o,z and -Bo.aj respectively, and the velocity field as recon- 
structed from this information (Fig. 3d). The comparison with UDV measurements 
has shown a good coincidence of the resulting velocity fields (see [55]). 

While GIFT is able to infer full 3D velocity fields by applying subsequently two 
different external magnetic fields, for the following case related to thin slab casting it 
can be reduced to the determination of the velocity components parallel to the wide 
faces of the mould 51 52 53J. Figure 4a shows the set-up of the corresponding model 
experiment Mini-LIMMCAST [S3] in which a liquid metal (in our case GalnSn) is 
poured from a tundish through a submerged entry nozzle into the mould. For this 
configuration it is sufficient to apply only one magnetic field by a single coil (see Fig. 
4a). The interaction of the fiow with the applied field produces induced magnetic 
fields that we measure at 2x7 positions at the narrow faces of the mould in order to 
reconstruct from them the velocity field. The mathematics of this inversion relies again 
on the minimization of the mean squared deviation of the measured magnetic fields 
from the fields according to the integral equation system (|18m9p . This minimization 
is done by solving the normal equations, whereby we use various auxiliary functionals 
which serve to ensure the divergence-free condition of the velocity, to enforce its two- 
dimensionality, and to minimize in parallel its mean quadratic value, weighted by 
some regularization parameter. 

Figure 4b and 4c show two examples of the reconstructed velocity fields resulting 
from this inversion. Note that due to the bubbling of argon into the liquid metal at 
the entrance of the submerged entry nozzle, quite different fiow structures can result, 
such as the desired double vortex flow (Fig. 4b), but also the unwanted single vortex 
flow as shown in Fig. 4c. 




Fig. 4. Application of GIFT to a flow problem related to the continuous casting of steel, 
(a) Experimental set-up: A liquid metal (GalnSn) is poured through a submerged entry 
nozzle (SEN) into the mould. Added argon gas leads to a highly chaotic flow in the mould. 
A magnetic held of about 2 mT is applied by a coil. The flow induced magnetic flelds are 
measured with 2x7 fluxgate sensors at the narrow faces of the mold. (b,c) GIFT reconstructed 
velocity field in the mould at two different instants, with a double vortex flow (b), and a 
single vortex flow (c). 
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4.2 Inverse problems at large Rm 

When extending the range of inverse problems from low to large Rm, we have to 
cope with highly nonlinear inverse problem. Some progress in its treatment has been 
made for restricted set-ups, for example by applying the so-called frozen-flux ap- 
proximation for the Earth's core which allows to determine (still with some appro- 
priate regularization) solutions for the tangential velocity field components at the 
core-mantle boundary from the measured time dependence of the radial magnetic 
field component [55^. Nowadays, such inversions are augmented by combinations with 
three-dimensional geodynamo modeling using methods like stochastic inversion [56] 
or data assimilation [5 7) . 

A complementary sort of restricted models is concerned with the determination of 
the radial dependence of the helical turbulence parameter a{r) from spectral proper- 
ties, under the (highly artificial) assumption that a is spherically symmetric. Solving 
such inverse spectral dynamo problems by means of an evolutionary strategy it was 
possible to obtain such a(r) profiles that lead to oscillatory dynamo solutions [55] . 
This model was later used to develop simple models of geomagnetic reversals that 
are just based on noise triggered jumps and relaxation oscillations in the vicinity of 
spectral exceptional points |59I60I61] . On this basis, a first attempt was further made 
to infer some of the most essential parameters of the geodynamo, such as its over- 
criticality and the turbulent resistivity from the temporal behavior and the statistical 
properties of field reversals \62\ . 




We have to admit, however, that apart from those special solutions with reduced 
dimensionality, inverse dynamo theory is still in its infancy. 



5 Conclusions and prospects 

In this paper we have discussed two ways of treating induction problems in magneto- 
hydrodynamics, one based on the differential equation approach, the other one based 
on the integral equation approach. We have applied both methods to various forward 
and inverse problems at small and large magnetic Reynolds numbers. 

Much work remains to be done in the field. Concerning the forward problem, 
the inclusion of time-dependent velocity fields can lead to surprising effects on the 
thresholds of self-excitation, based on non- normal growth (63] or parametric resonance 



As for the inverse problems at low Rm, externally applied fields with different 
frequencies could in principle be employed for better resolving the velocity field in the 
depth of the fiuid. Likewise, the use of AC fields might give a chance to combine CIFT 
with other tomographic methods, such as the Mutual Inductance Tomography for the 
identification of conductivity gradients, as they are connected with the presence of 
gas bubbles in liquid metals [55] . 

Inverse dynamo problems will remain a tough nut in the future, with some inter- 
esting prospects for inverse spectral methods and data assimilation. 
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